A <- read.table(file='cs-umass-matrix.dat');
A <- as.matrix(A)

ind=which(rowSums(A)>0)
A = A[ind,ind]
n=nrow(A)

#save original indices to ind
ind=c(375134,22484,104308,145622,370659,40421,206948,95996,219784,302146,186246,395649,34697,340131,135653,102074,71189,147703,128142,394741,45071,33796,330503,161514,709187,222416,140084,219592,137179,41394,40717,95006,102507,46426,46686,134726,131039,159910,39551,209598,317365,37055,12087,40261,78272,78116)
